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ABSTRACT 

We present a detailed model atmosphere analysis of high-dispersion and high signal-to- 
noise ratio spectra of the heavily polluted DAZ white dwarf GALEX J1931+0117. The 
spectra obtained with the VLT-Kueyen/UV- Visual Echelle Spectrograph show several 
well-resolved Sin spectral lines enabling a study of pressure effects on line profiles. We 
observed large Stark shifts in silicon lines in agreement with theoretical predictions 
and laboratory measurements. Taking into account Stark shifts in the calculation of 
synthetic spectra we reduced the scatter in individual line radial velocity measurements 
from ~ 3 to < 1 kms -1 . We present revised abundances of O, Mg, Si, Ca, and Fe 
based on a critical review of line broadening parameters and oscillator strengths. The 
new measurements are generally in agreement with our previous analysis with the 
exception of magnesium with a revised abundance a factor of two lower than previously 
estimated. The magnesium, silicon and iron abundances exceed solar abundances, but 
the oxygen and calcium abundances are below solar. Also, we compared the observed 
line profiles to synthetic spectra computed with variable accretion rates and vertical 
abundance distributions assuming diffusive steady-state. The inferred accretion rates 
vary from M = 2 x 10 6 for calcium to 2 x 10 9 gs -1 for oxygen. We find that the 
accretion flow must be oxygen-rich while being deficient in calcium relative to solar 
abundances. The lack of radial velocity variations between two measurement epochs 
suggests that GALEX J1931+0117 is probably not in a close binary and that the 
source of the accreted material resides in a debris disc. 
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1 INTRODUCTION 

GALEX J193156.8+011745 (GALEX J1931+0117, there- 
after) is a hydrogen-rich whit e dwarf discovered by 
IVennes. Kawka. fc Nemethl (|2010al ) and characterized by an 
opulent heavy-element line spectrum and an infrared ex- 
cess. The original low-resolution spectrum obtained with 
the New Technology Telescope (NTT) at La Silla Observa- 
tory showed a strong MgnA4481 doublet and weaker silicon 
lines. Follow-up echelle spectroscopy obtained with the Very 
Large Telescope (VLT)-Kueyen enabled a detailed abun- 
dance study. The near-solar abundances of oxygen, magne- 
sium, silicon, calcium and iron bear the signature of an exter- 
nal supply of material accreting onto the surface of the white 
dwarf. Based on available data, Vennes et al. concluded that 



* Based on observations made with ESO telescopes at the La 
Silla Paranal Observatory under programme 283.D-5060. 
f E-mail: vennes@sunstel.asu.cas.cz (SV); kawka@sunstel.asu. 
cas.cz (AK); nemeth@sunstel.asu.cas.cz (NM) 



the supply may originate from a close, sub-stellar companion 
or from a cool debris disc. 

The presence of heavy elements in hydrogen-rich 
white dwarfs has variously been interpreted as intrin- 
sic to the white dwarf, or as extrinsic, i.e., supplied 
by th e interstellar medium (|Dupuis. Fontaine, fc Wesemaell 
Il993t ). by a ne arby compan i on as in pos t -com mon en- 
velope systems dDebesl l200fj; iKawka et al 
a debris disc IjZu ckcrman et al 



20081). 



20031 : iKilic et al l [200 



>oi 



iFarihi. Zuckermam fc Becklinl 120081 ). However, accretion 
from the inte rstellar medium is unlikely because of sup- 
ply shortages (|Farihi et al.ll2010al ). In the extrinsic scenar- 
ios, the elements are accreted and d iffused in the atmo- 
sphere and envelop e of the star (see iFontaine fc Michaudl 
1979; Kocstcr 2009). An intrinsic, or internal, reservoir of 
heavy elements is also possible, but in either scenario a 
self-consistent solution of the diffusion equation must ex- 
plore the effect of radiative accelerat i on on trace elemen ts 
(|Chaver. Fontaine, fc Wesemael|[l995l : IChaver et al.lll995l ). 
As a class, the polluted DA white dwarfs, or DAZs, 
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are often defined by the detection of the Call HfcK dou- 
blet i n optical spectra of cooler objects (|Zuckerman et all 
2003; iKoester et al.l 120051 ). or by the detection of car- 
bon a nd silicon in t he ultr aviolet spectra of warmer ob- 
jects (|Dupuis et all l2009al lbh. Exceptionally, MgllA4481 
is, so far, only detect ed in a handful of warm whit e 
dwarfs such as EG 102 (jHolberg. Barstow. fc Greenlll997l ). 
GALEX J1931+0117, and two warm white dwarfs from 
the Sloan Digital Sky Survey (S PSS) that show evi- 
dence of dusty and gaseous discs |G~ansicke et al.l 120061 ; 
iGansicke. Marsh, fc Southworthl \200H ). The presence of a 
large concentration of magnesium in the last two objects 
helped establish a strong link between heavy element pol- 
lution and dense circumstellar environments. Moreover, an 
infrared excess that cannot otherwise be explained by a cool 
companion, may be attribu ted to a dust ring as in the case 
of the white dwarf G29-38 (|Graham et al.lll990t ). 

We present new high-dispersion spectroscopic observa- 
tions that help elucidate the nature of the peculiar abun- 
dance pattern in GALEX J1931+0117 in support of the 
debris disc model. Our observations are presented in Sec- 
tion 2, and a detailed line profile and abundance analysis is 
presented in Section 3. In Section 3.1 we compile and evalu- 
ate line broadening parameters and oscillator strengths. In 
Section 3.2 we describe steady-state diffusion models and 
predicted abundance gradients in the stable atmospheres of 
hot white dwarfs. In Section 3.3 we present revised abun- 
dances based on a review of atomic data, we assess non-local 
thermodynamic equilibrium (non-LTE) effects on measured 
abundances, and we constrain the accretion rates on the at- 
mosphere of GALEX J 193 1+0 11 7. Finally, we conclude in 
Section 4. 



Table 1. Additional UVES line identifications. 
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a Continued from lVennes. Kawka. fc Nemethl d2010al l . 
b Blended with Sill A5056.317. 

DA white dwarf is not in a close binary system. The average 
of the two measurements is vn a = 37.2 ± 0.6 kms -1 . 

Tabic [T] u pdates the spectral line iden tifications pre- 
sented in lVennes. Kawka. fc Nemethl l|2010al 'l. A set of nine 
Fell lines were observed allowing us to improve our ear- 
lier iron abundance measurement. Two strong Sin lines are 
added to our list, but their line velocities deviate signifi- 
cantly from the rest velocity, set by Ha, by 6 and 7 kms -1 . 
A closer examination of other Sill lines observed at epoch 
1 also reveals similar radial velocity shifts that shall be in- 
vestigated in Section 3.3.1. The sulfur lines SnA5014.069, 
A5032.447, and A5432.815A were not detected. 



2 OBSERVATIONS 

Following the identification of GALEX J1931+0117 as a 
DAZ white dwarf we obtained a series of echelle spectra 
using standard settings with the UV- Visual Echelle Spec- 
trograph (UVES) attached to the VLT-Kueyen. All spec- 
tra were obtained with the slit width set at 1", offering 
a resolving power ii w 46000. The exposure times were 
1450 s. The first series of spectra were obtained on UT 2009 
Nov 12 (epoch 1, thereafter) with the dichroic #2 and the 
"HER_5" and "BK7_5" filters on the blue and red arms, re- 
spectively. The blue spectrum was centred at 4370A and 
covered the range 3757-4985A, and the red spectra were 
centred at 7600A and covered the ranges 5698-7532 and 
7660-9464A. The first series of spec tra were analysed by 
us (IVennes. Kawka. fc Ncmcth 2010a). The second series of 
spectra were obtained on UT 2010 Mar 15 (epoch 2, there- 
after) with the dichroic #1 and the "HER_5" and "SHP700" 
filters on the blue and red arms, respectively. The blue spec- 
trum was centred at 3900A and covered the range 3290- 
4518A, and the red spectra were centred at 5800A and cov- 
ered the ranges 4788-5750 and 5839-6808A. 

We remeasured the Ha radial velocity by fitting Gaus- 
sian profiles to the well-exposed, narrow line cores obtained 
at epochs 1 and 2, i>h q ,i = 37.6 ± 0.8 kms -1 and vua,2 = 
36.8 ± 1.0 kms -1 , respectively. The two measurements are 
identical within errors. A lack of variations suggests that the 



3 ANALYSIS AND DISCUSSION 

We adopted the mod el atmo sphere analysis of 
IVennes. Kawka. fc Nemethl (|2010al ) that was based 
on non-LTE calculat i ons made with tlu sty-synspec 
jHubenv fc Land 1 19951 ; lLanz fc Hubenvlll995l ). All models 
were computed at T eS = 20890 K and logp = 7.90. 
I Vennes. Kawka. fc Nemethl l|2010al i found the difference 
between the best-fitting (T e ft, log g) to the Balmer line 
series using non-LTE heavy-element blanketed model atmo- 
spheres and the best-fitting solution using pure-hydrogen 
LTE models to be within statistical errors (AT c ff = +120 
K, Alog<? = —.04). Heavy elements contribute less than 
0.1% of the total electron density and their combined 
opacities are relatively modest in this range of effective 
temperatures compared to the strong Lyman and Balmer 
lines and continua. 

The model atmosphere structures and synthetic spec- 
tra employed in the line profile analysis were computed as- 
suming non-LTE, and with varying abundances of oxygen, 
magnesium, silicon, calcium and iron. The spectral synthe- 
ses are computed with different atomic data sets described 
below (Section 3.1). Also, we computed spectral syntheses 
with depth-dependent trace element abundances that were 
prescribed by diffusion theory (Section 3.2). Finally, we de- 
termine the extent of non-LTE effects in the atmosphere of 
GALEX J1931+0117 by comparing our non-LTE analysis to 
LTE results (Section 3.3.2). 
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Table 2. Sill oscillator strengths (/ij). 



Table 3. Sill Stark widths (FWHM) and shifts. 
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a Compiled bylKurucz fc Belli lll995h . 
l Artru et"al] jl98lf). 



i\i li u L'l ai. t lyoi 
Matheron et all j200lh . 



3.1 Atomic data 

Our original spectral syntheses were com puted with the 
line li st available on the CD-ROM No. 23 of iKurucz fc Belli 
|l995l fl. For each ion, we now compare the oscillator 
strength s (/ij) and Sta r k line broadening parameters (T) 
listed bv lKurucz fc Belli l|l995h to the best available theoret- 
ical and experimental data. Data on line oscillator strengths 
are also available at the National Institute of Standards and 
Technology (NIST)Q 

We noted that the isotopic shift in the 

26- 24 MgllA2798 doublet is o nly +0.85 kms" 1 

ijDrullinger. Wineland. fc Bergquistl Il980t l and, there- 
fore, we do not expect observable effects on the MgllA4481 
doublet in 6 km s _1 -resolution spectra. Similarly, other 
isotopic shifts (e. g., 30 Si/ 28 Si) may be neglected in the 
present study (see lBerengut. Dzuba. fc Flamba um 2003). 



3.1.1 Sin 

The rich silicon line spectrum in GALEX J1931+0117 
prompted a detailed review of available data. Table [5] lists 
and compares /ij from popular data compilations (CD23 
and NIST) to homogeneous theoretical or experimental 
data sets. We noted considerable variations in Sill oscil- 
lator strengths, in particular in the A3862 triplet and in 
A5041.024. Discrepancies of the order of 40 to 50% would af- 
fect individual abundance measurements in equal measures. 
For example, the ratio of CD23 data to experimental /ij val- 
ues is 1.06 but varies with a standard deviation a =38% 
while the NIST data that are largely based on the exper- 
imental values vary by only 8 % with an average ratio of 
0.98, and the theoretical data (|Artru et al.lll98lf ) vary by 
23% with an average ratio of 0.98. The adoption of one data 
set over another will not have a large effect on the aver- 
age abundance, but individual line measurements are less 
reliable. 
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a Compiled bv lKurucz 
^ Lanz, Dimitriicvic, 



and T = 20 X 10 3 K. 



Bell l|l995l) . 
Artrul dl988h calculated at n e 



10 1 



At n e 



10 1 



l2002j;lLesaeell2009h . 



and T = 16 - 20 X 10 3 K l lGonzalez et all 



Strong saturated lines are sensitive to line broaden- 
ing parameters. Table [3] lists full-width at half-maximum 
(FWHM= 2w) Stark widths an d shifts due to electro n im- 
pacts for strong Si II optical lines. IKurucz fc Belli (|l995l ) tab- 
ulate the circular frequency per electron F that we converted 
into the FWHM at n e — 10 lr cm~ 3 using the formula: 



(1) 



1 CD23, accessed at http: //www. cia. harvard. edu/amp/ampdata/kuruc 

2 Accessed at http://www.nist.gov/physlab/data/asd.cfm 



1 A 2 
2w — — F — n e , 

2tt c ' 

where c is the speed of light. lLanz. Dimitriievic. fc Artrul 
|l988r i tabulates FWHM values for electron and proton im- 
pacts separately. For most lines the proton contribution to 
the total width is » 10% up to 15%. The FWHM from 
IKurucz fc Bell (|l995h are on average 43% larger than the 
experimental v alues with a standard deviation of 4 4% while 
the values from lLanz. Dimitriievic. fc Artrul l|l988h are only 
6% lower than the experimental values with a standard de- 
viation of 21%. 

The excellent ag reement between 

lLanz, Dimitriievic. fc Artrul (|l988l ) and the experiments 
prompted us to explore two options in computing detailed 
silicon line spectra. In option 1 we adopt ed the silicon 
oscillator strengths from lArtru et al.l (|l98ll ) and the line 
broadening parameters for electr ons and protons from 
lLanz. Dimitriievic. fc Artrul (|l988l h The line broadening 
parameters are tabulated at 5, 10, 20, and 40xl0 3 K with 
estimated uncertainties of less than 20% at n e — 10 lr cm -3 
(or a depth tr ~ 2 in the atmosphere). Also, the effect of 
Stark shifts ar e inclu ded using the experimental data of 
iGonzalez et~all (120021 ) although we have no information on 
the scaling of d e with temperature or on the magnitude 
of Stark shifts due to ions (protons in this case). The 
effect of Stark shifts is apparent in radial velocity measure- 
ments of the Si ii A A5041-5055 and AA5967-5978 multi plets 
(see Table [[] and IVennes. Kawka. fc Nemethl bOlOal ). In 
option 2, we employe d the data (/ij and F) compiled by 
famtt flaifctjliaii (I1995I ). and neglected the effect of Stark 
shifts. 
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Table 4. Mgn Stark widths (FWHM) and shifts. 



Table 5. Call Stark widths (FWHM) and shifts. 
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4481.126 


1.09 


2.50 -0.0883 


2.95 -0.085 


3933.663 


0.248 


0.300 -0.048 


0.286 -0.14 



a Estimated using the classical formula and n a g = 4. 

b IDimitriievic fc Sahal-Brechotl l|l995l l calculated at n e = 1.35 X 

10 17 cm" 3 and T = 16.9 X 10 3 K. 

c Measured at n e = 1.35 X 10 17 cm" 3 and T = 16.9 X 10 3 K 
llDienize. Sreckovic, fc Bukvio!l2005r , lLesage|[2009h . 

3.1.2 Mgll 

Unlike silicon, the magnesium abundance measurement is 
based on a single doublet. Th e red MgllA7989 multiplet 
shows evidence of Stark shifts l|Vennes. Kawka. fc Nemethl 
l2010al ) although the spectrum is particularly noisy in the 
vicinity of the multiplet. On the other hand, the strong 
MgnA4481 doublet is well exposed. Because of its strength, 
the line profile is particularly sensitive to broadening param- 
eters, although quoted oscillator stre ngth values are consis- 
tent (NIST and lKurucz fc BellHl995Tl within < 2%. 

The compilation of IKurucz fc Belli (|l995h does not pro- 
vide a Stark width for the MgllA4481 do ublet which is then 
estimated in SYNSPEC using the formula (|Castellill200"5h 



r = icr 8 n c 5 ff , (2) 

where F is expressed in rad s _1 cm -3 and n c g is the effective 
quantum number in the hydrogenic approximation 

2 13.595 



2 

"-off 



(z + iy 



E t 



(3) 



where Ei is the ionization energy (eV) of the upper level of 
the transition and Z is the ionic charge of the trace element. 
The F value is then converted into the FWHM 2w. Table [J 
shows that theoretical and experimental values of the widths 
and shifts due to electrons are in agreement but that the es- 
timated width using the classical formula is in error. The 
Star k shift due to electrons compares well with experiments 
(see lLesagell2009l ). although the shift due to protons is pre- 
dicted to d ominate llDimitriievic fc Sahal-Brechotlll995l ). In- 
cidentally, IDimitriievic fc Sahal-Brechotl 1 19951 ) also predict 
a pressure shift in the Mg IIA7989 multiplet of opposite sign 
and a factor of ~8 larger than in the MgllA4481 doublet, in 
qu alitative agreement with ob s ervatio ns reported here and 
by I Vennes. Kawka. fc Nemethl (|2010al ). 

In step with our investigation of silicon lines, we also 
developed two options for the analysis of the MgnA4481 
doublet. In option 1 we adopted the NIST oscillator 
strengths and the line broa dening and shift parameters 
(proto n and electron) from IDimitriievic fc Sahal-Brechotl 
Jl995l). In option 2, w e adopted the oscillator strengths from 
IKurucz fc Bell (|l995l ) and the classical line broadening pa- 
rameter while neglecting pressure shifts. 



3.1.3 Call 

Table compares theoretical and experimental results 
for the Ca IIA3933 resona nce line. Star k widt hs from 
various sources including IKurucz fc BelH (|l995h are m 



Compiled bv lKurucz fc Belli l|l995h and accessed at http: 
//www. cf a. harvard . edu/amp/ampdata/kurucz23/sekur .html. 
6 IDimitriievic fc Sahal-Brechotl dl992h calculated at n e = 1.76 X 
10 17 cm" 3 and T = 43 X 10 3 K. 

c At n e = 1.76 x IP* 7 cm" 3 and T = 43.0 x 10 3 K 
dSreckovic fc Dienizejri993h . 



agreement although the pr e ssure shifts 
IDimitriievic fc Sahal-Brechotl ll 19921) 



calculated by 
l f actor of 3 

smaller than measured (|Sreckovic fc Dienizd [l993). Be- 
cause the measured CallA3933 sh ift is relatively small (see 
IVennes. Kawka. fc N emcth 2010a) it may not be measurable 
in our spectra and the lack of accuracy in these parameters 
would not affect the abundance and radial velocity analysis. 

As part of option 1 we adopt the theoretical broad- 
ening and shift parameters of IDimitriievic fc Sahal-Brechotl 
ll 19921) and in option 2 we adopt again the compilation of 
IKurucz fc "Bell (l995l ) without pressure shifts. 



3.1.4 Ol and Fell 

The analysis of weak spectral lines is sensitive to the as- 
sumed oscillator strengths but it is less sensitive to broaden- 
ing parameters, and we will not distinguish the two options. 

The O IA 77 73 o scillator strengths tabulated by 
IKurucz fc Bell l| 19951 ) and NIST differ by wl0% 
which implies a minimum uncertainty in oxy- 
gen abundance measurements of the same order. 



iBen Nessib. Lakhdar. fc Sahal-Brechot d 19961 ) compared 
calculated Ol Stark width to iGriem (| 19741 ) and to ex- 
periments. They found excellent agreement between their 
calculations and the experimental values although they do 
not provide data for the O IA7 773 multiplet. The Stark 
width tabulated bv lGrieml l| 19741 ) and lKurucz fc Belli (|l995l ) 
agree to within a few percent at T = 20000 appropriate for 
our analysis. We ad opted the oscilla t or str engths and Stark 
widths tabulated bv lKurucz fc Belli (|l995l ). 

Five of the Fen lines detected in GALEX J1931+0117 
are listed at NI S T wh ile atomic data for all lines are listed by 
IKurucz fc Bell (|l995l ). In this case, we adop ted the oscilla- 
tor str engths and Stark widths tabulated bv lKurucz fc Belli 
|l995l ). 



3.2 Diffusion in white dwarf atmospheres 

A tenet of the accretion-diffusion scenario is that photo- 
spheric abundances reach a steady-state level in the line- 
forming region. The abundance gradient is assumed to 
be zero in convective layers, but it is a function of the 
diffusion time scale in radiative layers as in the case of 
GALEX J1931+0117. 

Assuming accretion steady-state the mass-continuity 
equation may be written as 
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C2 



^=M 2 



1 n,2 



4nR 2 A 2 m 



(4) 



where c 2 is the abundance of element "2" relative to the 
main constituent "1", M 2 is the mass accretion rate of el- 
ement "2" (gs -1 ), R is the stellar radius (cm), A 2 is the 
atomic weight of element "2" (a.m.u), 



Tl,2 



m/{pvi l2 ) 



(5) 



is the diffusion time scale of element "2" (s), m and p are 
the mass loading (gem -2 ) and density (gem -3 ) in the at- 
mosphere, and t>i,2 is the diffusion velocity (cms -1 ) given by 
the diffusion equation. Assuming that element "1" is fully 
ionized hydrogen (Ai = f, Z\ = 1) and that c 2 « 1, then 
the diffusion equation fsee lVennes et al"1ll988T l simplifies to 



«1,2 = -Dl 



V c 2 dr 



A 2 



1 dc 2 

Ox 

is the diffusion 
relative to the 



m p9 . , m v9 
+ A 2 



(6) 



kT kT 

coefficient of trace el- 
main constituent (see 



where Di )2 
ement "2" 

IVennes. Kawka. fc NemethllioiQbl ) 

The first term on the right-hand side is the actual abun- 
dance gradient. We show in Section 3.3.3 that c 2 is a slowly 
varying function of density, c 2 oc p a and \a\ < 0.5, therefore 



(7) 



fdc2 < 1 dp dp 
c 2 dr ~ p dr dm' 

where dm — —pdr. The second term on the right-hand side 
of the diffusion equation is proportional to the surface grav- 
ity g and the last term is proportional to the radiative accel- 
eration g r . Employing the run of density p and temperature 
T as a function of mass loading m in a model atmosphere 
appropriate for GALEX Jf 93f+0f f 7, it can be shown that 
gravitational settling acts unimpeded as long as c 2 << 1, 
i.e., that 



dp 




dm 


« 





(l-A 2 ) 



m p g 



kT 



(8) 



The abundance gradient term is < 2% of the gravity term 
and does not effectively co unteract the gravitational pull. 
Simil a rly, radiative forces jChaver. Fontaine, fc Wesemaell 
1 19951 ; IChaver et al.l Il995l ) are relatively weak at T e g ~ 
20 x 10 3 K and cannot raise the abundances to levels observed 
in GALEX J1931+0117. The silicon abundance potentially 
supported by radiation is at least 3.5 orders of magnitude 
below the observed abundance, and the radiative forces on 
oxygen, magnesium, calcium, and iron are negligible. 

In summary, unless a diffusion steady-state cannot be 
reached, the vertical abundance distribution in the atmo- 
sphere of an accreting white dwarf follows a relation of the 
form C2 oc Ti l2 /m and may depart significantly from a homo- 
geneous distribution. Longitudinal homogeneity is assumed 
although abundance spots may in principle occur. 



3.3 Line profile analysis 

Detailed line profiles are computed and convolved to the 
instrumental resolution (R « 46000). The separate routines 
were introduced into the spectral synthesis code to treat 
Stark shifts self-consistently. The Voigt profiles, that already 
include a treatment of Doppler width and damping functions 
(r nat + Tstark) have been modified to include the effect of 
Stark shifts. The line centre is shifted linearly with electronic 
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Figure 1. Abundances measured in LTE with individual lines 
and using two separate atomic da ta sets. Set number 1 em- 
ploys line broadening calculations of lLanz. DimitriieviCj fc Artrul 
l|l988h and oscillator strengths from lArtruet*lrt* I \l98Tf) and in- 
cludes the ef fect of Stark shifts. Set number 2 employs data 
compiled by iKurucz fc Belli l ll995l) but excludes Stark Shifts. 
Twenty abundance measurements (full circles) with set num- 
ber 1 are shown in the upper panel, while the same abundance 
measurements using set number 2 are shown in the lower panel 
al ong with the non-LTE mea s uremen ts (open circles) collected 
bv lVennes. Kawka. fc Nemethl (2010ah and showing the extent of 
non-LTE effects on the silicon abundance measurements. Solar 



abundance (log n(Si)/n(H) 



-4.4 9) is shown with dotted lines 



llAsplund, Grevesse. fc Sauva and the measured averages 

are shown with dashed lines. Standard deviations (Her) for each 
data set are shown in the lower right corner (see text). 



and proton density and proportionally to the Stark shift d 
tabulated at a reference density. Table [6] lists the revised 
abundances. 



3. 3. 1 Line broadening and pressure shift 

In the present analysis we assume both longitudinal and 
vertical abundance uniformity. The lack of abundance vari- 
ations between two observation epochs favours a uniform 
spread of trace elements across the surface although addi- 
tional observations are required to rule out a mere coinci- 
dence in rotational phases. Figure [1] shows abundance mea- 
surements obtained for nine Si II spectral lines in 20 separate 
measurements and Figure [2] shows sample line profile fits. 

The revised silicon abundance is —0.05 dex (or 10%) 
lower than estimated bv lVennes. Kawka. fc Nemethl ([gOlOa) 
mainly because of the additional lines included in the anal- 
ysis. In fact, the abundance measurements obtained assum- 
ing option 1 were on average 7% larger than with option 
2. However, as proposed in option 1, adopting the silicon 
oscillator strengths and line-broa deni ng paramete r s from 
lLanz. Dimitriievic. fc Artrul (|l988h and lArtru et all (|l98lf ). 
respectively, reduced the scatter in abundance measure- 
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Figure 2. Line profile fits of spectra obtained at epoch 2 and cor- 
responding to abundances depicted in the top panel of Figure [l] 



merits based on 20 separate measurements. The scatter 
in abundance measurements obtained with option 1 was 
a — 0.08 while the scatter obtained with option 2 was 
larger, a — 0.13, while the measurements showed a definite 
wavelength-dependent trend (Fig. [TJ . 

Moreover, the inclusion of Stark shifts in line profile 
modelling reduced the scatter in Sill line velocity mea- 
surements. The standard deviation of the measurements is 
a v — 2.8 without shifts down to a v = 0.9 kms -1 with Stark 
shifts included. The average line velocity is v = —0.8 relative 
to tiHo = 37.2 kms -1 without Stark shifts and 0.0 kms -1 
with Stark shifts included. The modelling of the spectral 
lines with the largest expected Stark shifts has been consid- 
erably improved by taking this effect into account (Fig. 0. 
Figure[3]compares Sill line profiles computed with and with- 
out Stark shifts to the observed spectra. The zero velocity 
is set by the Ha line core (+37.2 kms -1 ). The line asymme- 
try in A5055 is due in part to the presence of a weaker Si n 
line over the extended red wing, but the asymmetric pro- 
file in A5041 is entirely due to the Stark shift. The velocity 
discrepancy notable for these two lines is entirely resolved 
using Stark shifted profiles. 

The calcium abundance was revised slightly down by 
—0.06 dex (or 14%) mainly because of a small increase in 
the adopted Stark width. The effect of pressure shift was 
not noticeable. We estimated the error on the abundance 
measurements by assuming an uncertainty of 25% on the 
broadening parameter. The corresponding uncertainty on 
the abundance is 0.05 dex or 10%. 

Our revised magnesium abu ndance is significantly lower 
than o riginally estimated by IVennes. Kawka. fc Nemetbl 
l|2010al ) entirely because of an upward revision of the Stark 
width by a factor of two in better agreement with experimen- 
tal measurements (Table 0}. We also estimated the abun- 
dance error assuming a 25% uncertainty, i.e., at least the 
difference between theory and experiments, on the broaden- 



Figure 3. Spectra obtained at epoch 2 are fitted with model line 
profiles including Stark shifts (full lines) or excluding the effect 
(dashed lines). The velocity shift in the observed line profiles is 
well reproduced by models including Stark shifts. 



ing parameter. The corresponding uncertainty on the abun- 
dance is 0.07 dex or 15%. 

The iron abundance, now based on nine separate line 
measurements, is, within errors, identical to the original 
estimate. The error was set to the standard deviation of 
the set of measurements. The procedure adequately propa- 
gates uncertainties in broadening parameters and oscillator 
strengths. However, variations of 25% on the Stark width 
had no effect on the line equivalent widths. Similarly, the 
oxygen abundance is identical to the original estimate, and 
the error was set to the corresponding uncertainty in oscil- 
lator strengths (10%). Variations in the Stark width had no 
effect on the equivalent widths because weak lines such as 
Fe n and O I lie on the linear part of the curve-of-growth 
and are insensitive to broadening parameters unlike strong 
Mgn, Sill, or Can lines. 

3.3.2 Non-LTE effects on abundance measurements 

We measured the abundances using the grid of LTE spectra 
and compared them to the non-LTE results. We found that 
non-LTE effects are the largest in highly-excited Sill lines. 
The average non-LTE shift on abundance measurements is 



A log Si/H = log Si/Hj, 



log Si/H LTE = +0.08 ± 0.03 



The abundance shift is w +0.05 for spectral lines with the 
lower levels 3s3p 2 ( 2 D) or 3s 2 4s( 2 S), and « +0.11 for highly- 
excited lines with the lower levels 3s 2 3d( 2 D) or 3s 2 4p( 2 P°). 

A similar comparison of magnesium, calcium, and iron 
abundances show that non-LTE effects are small: 

AlogMg/H = -0.02 

AlogCa/H = -0.02 

AlogFe/H = +0.01 
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Finally, the abundance of oxygen measured with highly- 
excited multiplet OlA7773 is relatively large: 

AlogO/H = -0.13 

In summary, non-LTE effects are notable in highly-excited 
lines of non-dominant species such as O I, but are negligible 
in strong lines of dominant species such as Can or Mgn. 
Highly-excited levels of Sill are also affected. 



3.3.3 Effect of abundance gradient on line profiles 

We investigated the effect of vertical abundance inhomo- 
geneities in the calculation of a spectral line syntheses. 
The depth-dependent abundance was calculated by fixing 
the accretion rate M and the stellar radius R and vary- 
ing the diffusion time-scale n,2 as a function of mass 
loading m as described in Section 3.2. We computed the 
time-scale and correspond i ng diff usion velocity following 
IVennes. Kawka. fe Nemethl l|2010al ). The effect of neutral 
species on element-averaged velocities was investigated by 
excluding them from the diffusion velocity calculations and 
comparing the resulting abundance profiles. 

Figure [4] (left) shows that in the line forming region 
the abundance ratio C2 = ri2/ni can be approximated by 
C2 oc p a where —0.44 < a < 0.54, or \a\ < 0.5. The simula- 
tions confirm that, in this context, the abundance gradients 
are not steep enough to counteract gravity (Section 3.2). The 
effect of neutral species on element-averaged diffusion veloc- 
ities are more marked in the cases of carbon and oxygen be- 
cause of the relatively high Cl and Ol ionization potentials. 
In the case of oxygen the assumed location of the photo- 
sphere bears heavily of the derived accretion rate and other 
elements may be affected as well. In order to remedy the sit- 
uation we computed new emerging line profiles parametrized 
by a constant accretion rate, i.e., an accretion steady-state, 
rather than a constant abundance throughout the photo- 
sphere. These non-LTE spectral syntheses were calculated 
by adopting the structure (T, p) of the best-fitting homo- 
geneous model and introducing the abundance distribution 
prescribed by Equation 4. The observed spectra were then 
fitted with the model grid and the accretion rates were mea- 
sured explicitly. 

Table [7| lists measured accretion rates for C, O, Mg, 
Si, Ca, and Fe. The quoted errors are statistical only and 
do not include systematic errors in modelling the diffusion 
flow. Calculated concentrations in the line forming region 
are inversely proportional to diffusion coefficients and errors 
would scale accordingly, i.e., an error of 50% on theoretical 
diffusion coefficients would propagate to calculated concen- 
trations. Figure [4] (left) shows the corresponding abundance 
profiles. The inferred abundances at three different locations 
in the photospheres (marked with arrows) are in good agree- 
ment with the abundance derived from homogeneous models 
in the case silicon, and they deviate slightly in the case of 
magnesium and iron. Figure [4] (right) shows model spec- 
tra compared to the observed spectrum at three accretion 
rates. Although abundance inhomogeneities may affect line 
profiles, in particular in the relative strength of line core 
and line wings, a comparison of ultraviolet to optical line 
strengths may be more revealing, particularly for ultraviolet 
oxygen lines. 

The inferred chemical composition of the accretion flow 



Table 6. Revised non-LTE abundances in GALEX J1931+0117. 



X 



n(X)/n(H) 



[X/H] a 





< 7 x 10~ 5 




< -0.54 


o 


2.4 ±0.3 x 10" 


-4 


-0.28 


Mg 


3.8 ± 0.6 x 10- 


-5 


0.05 


Si 


5.8 ± 1.0 x 10- 


-5 


0.25 


Ca 


7.7 ±0.8 x 10" 


-7 


-0.42 


Fe 


3.7 ±0.8 x 10" 


-5 


0.12 



a [X/H ]=Iogn(X)/n(H) 
^ From Venncs, Kawka, & 



-Iogn(X ) /n(H) 
Nemethl feOlOafi 



Table 7. Accretion rates for selected elements. 



X 


log M (gs" 1 ) 


[X/Si]^, 


c 


< 8.0 


-0.8 





9.30 ±0.05 


0.1 


Mg 


7.75 ±0.06 


-0.5 


Si 


8.25 ±0.07 




Ca 


6.39 ±0.04 


-0.9 


Fe 


8.45 ±0.09 


0.0 



a [X/Si]= logn(X)/n(Si) - logra(X)/n(Si) 

may be estimated by calculating the abundances relative to 
silicon 



n(X) 



n(Si) 



M(X) / M(Si) 
A(X) \ A(Si) 



(9) 



where X refers to C, O, Mg, Ca, and Fe. Table [7] list the 
corresponding abundance ratios relative to solar ratios. The 
results show that oxygen must be supplied in larger propor- 
tions that afforded by a solar composition, and that silicon is 
also supplied in larger proportions than carbon, magnesium, 
or calcium. Iron is apparently supplied in solar proportions 
relative to oxygen. 

In summary, a line profile analysis parametrized by ac- 
cretion rates rather than constant abundances reveals an 
oxygen enrichment in the accretion flow. Diffusion calcula- 
tions show the dominant effect of neutral species on abun- 
dance profiles. 



4 DISCUSSION AND CONCLUSIONS 

We presented a model atmosphere analysis of the high- 
metallicity white dwarf GALEX J1931±0117. Based on 
a critical review of available atomic data, their uncer- 
tainties, and their impact on the accuracy of abundance 
measurements, we revised the ab undances presented by 
IVennes. Kawka. fe Nemethl (|2010al ). The variations are not 
significant with the exception of magnesium which is re- 
vised downward by a factor of two. The abundance pattern 
shows that the abundance of magnesium, silicon, and iron is 
near or above the solar abundance, while the abundance of 
carbon, oxygen and calcium is below the solar abundance. 
A line profile analysis performed using vertical abundance 
distributions obtained by solving the steady-state diffusion 
equation shows that the accretion flow is rich in oxygen and 
silicon while it is depleted in other elements. The calcula- 
tions show that although the oxygen abundance is below so- 
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log p (g cm- 3 ) v (km s -i) 

Figure 4. (Left) vertical abundance distribution of elements in the atmosphere of GALEX J1931+0117 computed assuming diffusive 
steady-state (Section 3.2) and typical accretion rates M in gs _1 (full lines). The effect of neutral species on element- averaged diffusion 
velocities and vertical abundance distributions is shown by excluding neutral species from the calculations (dashed lines). The profiles 
are compared to measured abundances from Table l6l (dotted lines). The locations of optical depth 2/3 at A = 1260A (rightmost) and in 
the line wings and center (leftmost) are shown with arrows. (Right) spectrum obtained at epoch 2 (jagged line) compared to model line 
profiles assuming Mg accretion rates M = 3.2, 5.6, and 10 X 10 7 gs _1 . 



lar in the line forming region, it must be supplied in larger 
quantity because of its short diffusion time-scale relative to 
other elements. 

On the other hand, we show that introducing Stark 
shifts in line profile modelling helped reduce the scatter in 
line velocity measurements to well within the expected ac- 
curacy of the echelle spectra. Investigations of the effect of 
pressure shifts on radial velocity measurements ar e relatively 
rare but their effect have been pre dicted (e.g., iHammondl 
ll989l ; lKrslianin fc Dimitriievic1ll992T ). A lack of radial veloc- 
ity variations between two epochs that are 123 days apart 
also rules out a close binary scenario for the origin of the 
accreted material. We are left with the possibility that the 
infrared excess belongs to a warm disc of debris material 
that accretes onto the white dwarf surface. 

T he iron abundance in GALEX J19 31+0117 is nearly 
solar |Asplund. Grevesse. fc Sauva 1 2005) and it is compa- 
rable to the iron abundanc e measured in much hotter white 
dwarfs (|Vennes et al.ll2006l ) . The resurgence of iron in cooler 
white dwarfs is also evidence that an external reservoir 
reached the surface of these objects. 

Simulated infrared spectral energy distributions for 
GALEXJ 1931+01 17 including a spectral synthesis of the 
white dwarf added to a warm disc show that the ob- 
served distribution jVennes. Kawka. fc Nemethll201o"ah can 
be matched by a warm disc (900 K) or a hot disc (1650 
K) near sublimation temperature. Necessarily, the inferred 
dus t temperature must be below sublimation temperatures 
(see[Loddcrs 200^) for the main constituents such as calcium 
(1659 K), or silicon and magnesium (1354 K). Therefore, the 
temperature as well as the size and state (gas or solid) of the 
disc is the result of an adjustment with the host white dwarf 



temperature and luminosity. The corresponding disc sizes 
are 90Qwd for the hot disc, and, compensating for its re- 
duced emissivity per unit area, 2200fiwD for the warm disc, 
where Qwd is the solid angle subtended by the white dwarf 
itself. The observed spectral energy distribution extends to 
the 2MASS measurements in the infrared and measurements 
at longer wavelengths are required to further constrain the 
model. 



A comparison with other warm DAZ white dwarfs sup- 
ports these conjectures. The DAZ PG 1015+161 (T eB = 
19300 K, logy = 7.9, logCa/H = -5.9. iKoester fc Wilkenl 
|2006| ) is similar to GALEX1931+0117 (T eS = 20890 K, 
loggr = 7.9, l ogCa/H = —6.1). Based on JHK and Spitzer 
observations. I Jura. Farihi. fc Zuckermanl |2007l ) inferred the 
presence of a disc with an inner temperature ranging from 
800 to 1000 K a nd emitting area from 700 to 12 00 Q. W d- 
In a similar case, iFarihi. Jura, fc Zuckermanl (|2009l ) inferred 
the presence of a disc with a blackbody temperature of 
1500K around PG1457-086 (T cff = 20400 K, \ogg = 8.0, 
logCa/H = —6.1). The great diversity in disc emissivity is 
illustrated by the case of HS0047+1903. Despite having a 
similar calcium abundance, no evidence of a disc is found in 
JHK and Spitzer observations |Farihi et al.ll2010bl ). 



In conclusion, our analysis of the abundance pattern in 
GALEXJ1931+0117 and the absence of a close companion, 
as well as a comparison with other, similar cases such as the 
DAZ PG1015+161 support the likely presence of a dusty 
disc around GALEXJ1931+0117. 
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